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C — , Abstract. We describe a general procedure for using number counts of any object to constrain 

the probability distribution of the primordial fluctuations, allowing for generic weak non- 
Gaussianity. We apply this procedure to use limits on the abundance of primordial black 
holes and dark matter ultracompact minihalos (UCMHs) to characterize the allowed statistics 
of primordial fluctuations on very small scales. We present constraints on the power spectrum 
. ^ and the amplitude of the skewness for two different families of non-Gaussian distributions, 

^ distinguished by the relative importance of higher moments. Although primordial black holes 

probe the smallest scales, ultracompact minihalos provide significantly stronger constraints 
on the power spectrum and so are more likely to eventually provide small-scale constraints 
on non-Gaussianity. 
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1 Introduction 

The power spectrum of primordial fluctuations is well measured on cosmological scales using 
the Cosmic Microwave Background (CMB) and Large Scale Structure [1—6]. These obser- 
vations provide compelling evidence that the fluctuations originated in an era of inflation. 
However, there is a great deal of new information waiting to be accessed in small-scale fluc- 
tuations and in higher-order statistics (non-Gaussianity). In this paper we explore both new 
regimes by examining how object number counts constrain the probability distribution of 
primordial fluctuations. By looking at very small objects (primordial black holes and ultra- 
compact minihalos), we constrain the power in fluctuations on smaller scales than the CMB 
and Large Scale Structure currently probe, k > 3 Mpc -1 . Since these objects are also very 
rare, we constrain the level of non-Gaussianity. 

Counts of rare objects are a useful probe of the primordial inhomogeneities in the grav- 
itational field and their evolution: different objects probe different scales, different cosmo- 
logical eras, and different particle physics and astrophysics. Primordial black holes (PBHs) 
have been used as a probe of small-scale power for many years [7]. Their abundance limits 
the allowed fluctuation power on very small scales, but the constraint is weak compared 
to the numbers expected from extrapolating a scale-invariant spectrum down from CMB 
scales [8, 9]. More recently, a much stronger constraint on small-scale power was obtained by 
considering ultracompact minihalos (UCMHs) of dark matter [10-12]. The only drawback of 
this approach is that the strongest limits [11] require dark matter to annihilate into Standard 
Model (SM) particles, which can then be sought using standard indirect dark matter detec- 
tion techniques [13]. Nonetheless, UCMH constraints on the small-scale power obtainable by 
microlensing [12] are still stronger than PBH constraints, and apply even for non-annihilating 
dark matter. 

The abundance of PBHs was also considered early on as a probe of non-Gaussianity [14] . 
In a non-Gaussian distribution, the number of rare objects is different from the Gaussian ex- 
pectation, and rarer objects are typically sensitive to higher moments of the distribution. 
Number counts pick up any deviation from Gaussianity in a model independent way and 
provide complementary constraints to those on the individual correlation functions, such as 



- 1 - 



the shape of the 3-point function in momentum space (the bispectrum). However, there are 
important limitations on what can be learned from number counts on both the analytic and 
observational sides. First, in general one only knows the non-Gaussian probability density 
function (PDF) approximately, usually in terms of a few of the lowest-order moments. This 
limits how far out onto its tail the distribution is known, and so limits the utility of looking 
at very rare objects to constrain particular models. Furthermore, the best controlled ap- 
proximations to weakly non-Gaussian PDFs are asymptotic expansions, so for a given level 
of non-Gaussianity one trusts the expansion only so far out onto the tail. The greater the 
level of non-Gaussianity, the more limited the range of utility of the asymptotic expansion. 
On the observational side, one must know the amplitude of fluctuations (the variance of the 
distribution) very precisely in order to find sharp constraints on the level of non-Gaussianity. 
The relationship between mass and the observable signature of the objects must be also be 
known accurately and precisely. Finally, extremely rare objects are not necessarily expected 
to be present in any finite sample (e.g. in any single survey or even our Universe), which 
makes drawing conclusions from their absence, or from very small number statistics, difficult. 
Nevertheless, constraints on the abundance of rare objects currently provides the only probe 
of the primordial fluctuations on very small scales, k ^ 10 4 Mpc -1 . On larger scales, 50 
Mpc -1 < k < 10 Mpc , //-type spectral distortions of the CMB also provide constraints 
[15-17]. We will return to the complementarity of these probes in the conclusions. 

Here we treat the analytic difficulties mentioned above carefully and develop a general 
prescription to use number counts to constrain weakly non-Gaussian primordial fluctuations. 
We apply this technique to derive bounds on non-Gaussianities with UCMHs and PBHs. 
Previous authors have used PBHs to constrain particular cases of the non-Gaussian PDF or 
particular inflation scenarios. Bullock and Primack [14], followed by Green and Liddle [18] 
and Ivanov [19], considered potentials with a localized feature that can both boost power 
dramatically on some scales and generate strong non-Gaussianity because of mode-coupling 
near the sharp feature. Pina Avelino [20] worked out constraints on a x 2 distribution for 
the fluctuations. Hidalgo [21] considered a non-Gaussian PDF from only the first term in 
the Edgeworth expansion (pi from Eq.(2.12)). Klimai and Bugaev [22] constrained two- 
field inflation models. Most recently Byrnes et al [23] considered constraints on weak or 
strong non-Gaussianity of the 'local' type [24], where the primordial curvature is lZ{x) = 
1Zg{x) + ^fNLR? G (x), where TZg(x) is a Gaussian random field. 

Our results generalize all previous work on weakly non-Gaussian models, demonstrating 
how constraints can be obtained from any object abundance and for any weakly non-Gaussian 
distribution. We use an extended Press-Schechter [25] approach where the abundance of 
an object is determined by the overdensity required for its formation (5 C = 5p c /p), the 
mass variance smoothed over the region that forms the object (<tr), and similarly smoothed 
higher-order moments of the density fluctuations. These variables can be straightforwardly 
computed for any object: one need only compute the collapse threshold and smoothed statis- 
tics (by propagating the primordial fluctuations forward to the appropriate time of formation 
using transfer functions and a growth factor). This method is a straightforward extension 
of previous work on an analytic, weakly non-Gaussian mass function for galaxy clusters 
[26]. However, here we provide more precise criteria for where the analytic expansions are 
reliable in terms of the rareness of the object (the threshold v c = 8 c /&r) and the level of non- 
Gaussianity. In addition, we allow a wider range of non-Gaussian scenarios by parameterizing 
deviations from Gaussianity with both the value of the dimensionless skewness (M3) and a 
rule for how the higher moments scale with respect to the skewness. Most previous work as- 
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sumes one particular choice, but physical mechanisms generating the primordial fluctuations 
do not all follow that pattern, and the distinction is important for interpreting constraints. 
By demonstrating the computational limits of the Edgeworth expansion, we also clarify why 
strongly non-Gaussian models must be evaluated on a case-by-case basis. 

Section 2 develops our formalism for describing weakly non-Gaussian distributions and 
determining the abundance of objects when the primordial fluctuations are described by a 
non-Gaussian distribution. Sections 3 and 4 apply this technique to PBHs and UCMHs 
respectively. We conclude by comparing current and future constraints on small scale power 
and non-Gaussianity on a range of scales from a variety of existing and proposed techniques. 

2 Calculating number counts using weakly non-Gaussian distributions 
2.1 Non-Gaussian Probability Distributions 

Assuming that objects form from regions of the density field above some threshold, we need 
a probability distribution P of perturbation amplitudes 5r and a threshold 5 C above which to 
count the objects. Here 5 Ct x gives the minimum density contrast required for a perturbation 
to later collapse into an object of type X G {PBHs, UCMHs}; <5 c ,pbh ~ 1/3, while <5 c ,ucmh 
is a scale-dependent quantity [11], although typically 0(1O~ 2 — 1CP 3 ) [11, 27]. Assuming a 
flat universe, the fraction of perturbations that will form an object of type X is 



where P is a function of quantities smoothed on the scale R associated with objects of mass 
M. For a generic PDF, the factor A can be determined by ensuring that counting all objects 
recovers the total mass density. Here we will just use the Press-Schechter factor of 2 (A = 1) 
as a good enough approximation since we are only concerned with a subset of objects 1 . In 
addition, we have approximated the integral as extending to infinity, even though for some 
rare objects (e.g. UCMHs) that are not quite as rare as others (e.g. PBHs) there will be 
a finite upper limit associated with 5 C of the 'next rarer' object. Extending the integral to 
infinity is a very good approximation as long as 5 C is dissimilar for different rare objects, as 
P is always a steeply-falling function of 5r. If the probability distribution is Gaussian, we 
have 



where v Ct x = 8 c> x I &R and or is calculated in the cosmological era relevant for formation of 
objects of type X. In general we will suppress the subscript X on /?, S c and i/ c , except where 
talking about a specific type of object. Our convention follows [11]: R refers to the comoving 
radius of the region that forms the object, k = 1/R is the corresponding wavenumber, and 
both S Ct x and <jr are computed when the region enters the Hubble horizon. This is an 
analytic technique, which should be adjusted to agree with numerical simulations before it 
is trusted to give precise predictions of object abundance. 

We can describe weakly non-Gaussian distributions as an expansion in moments about 
a Gaussian. A generic expression, generalizing the Edgeworth expansion, was written down 

1 This factor would be more important if one had a particular model for how the primordial statistics varied 
as a function of scale and combined constraints on that model from objects on several scales. 






(2.2) 
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by Petrov [28, 29] 



In this expression, the second summation is over all sets of integers {k m } that satisfy 

ki + 2k 2 + • • • + sk s = s ; (2.4) 

H n (u) are Hermite polynomials 2 H n (v) = {—l) n e v2 / 2 -^e~ v2 l 2 ; and r = k\ + k<i + • • • + k m . 
The "reduced cumulants" are defined as 

Sn,R = ,1^1-1 ' ( 2 - 5 ) 
\°R/c 

To decide how to organize the terms in this series and truncate it somewhere, we need 
additional input about the relative importance of higher order moments for the physical case 
under consideration. The choices are more naturally phrased in terms of the connected part 
of the dimensionless moments (the cumulants) of the real space density contrast, smoothed 
on scale R: 

^.* = (ffe- < 2 - 6 > 

The moments A4 n ,R can equally well be thought of as functions M n {k) of the wavenumber 
k corresponding to the smoothing scale R. We will use these two notations interchangeably, 
and sometimes we do not write out the scale-dependence explicitly at all. 
Two choices motivated by particle physics [30] are hierarchical scaling 



M h n oc \TP£\ , (2.7) 

and feeder scaling 

M l n oc l n , n > 3, (2.8) 

where V-r, is the amplitude of fluctuations in the primordial curvature, (7£ 2 (x)) = f ^-V-ji(k), 
and I is a parameter indicating the strength of the interaction that sourced the non-Gaussianity. 
In specific models one can work out the coefficients relating X oc /jvx, where /jyx is the pa- 
rameter for which WMAP quotes constraints on non-Gaussianity. Using simple local models 
with each scaling to determine representative combinatorics 3 , we can express the typical size 
of all higher order moments (n > 3) in terms of the third moment .M3 by 

Hierarchical M\ = n\ 2 n " 3 f ^ J (2.9) 

Feeder M { n = (n - 1)! 2 n ~ l (^f) " /3 . (2.10) 



2 Note that there are two conventions for naming the Hermite polynomials, and these are sometimes denoted 
He„. 

3 We use the local ansatz 1Z(x) — 1Zg(%) + § fNL{1ZG(x) 2 — {T^-g{x) 2 }] for the hierarchical scaling (in that 
case X — /jvl)- We use a two field extension 1Z(x) — <f>o + <?g + §/zvl[o"g(20 2 — ( a G{x) 2 )], with JnlPV 2 S> 1 
for the feeder scaling (then X — f'ff/ = /nlVct/Tk )■ 
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These are of course just two examples of how the moments could scale. For currently devel- 
oped inflation models, they are representative of the full range of well motivated possibilities. 
Some models (quasi-single field inflation [31]) are a hybrid between the two scalings, or have 
additional numerical coefficients (e.g. two-field local models). Any skewness in P(v) must be 
positive (assuming that higher moments are sub-dominant) for non-Gaussianities to increase 
the production of rare objects like PBHs, as the positive tail of the distribution must be 
enhanced relative to the core. 

Then Eq. (2.3) can be written in terms of the dimensionless smoothed moments J^A.n R 
and the series can be organized in a form appropriate for any scaling. Making use of Eq. 
(2.4), we find 



dv 



V/27T 



s=l{k m } m=l m VV ^ > / I 



(U ' e-» 2 / 2 (l+ Pl {v,R)+p 2 (v,R) + ...) 



2tt 

For the two scalings above, we have 

p[ h \v,R) = ~^H 3 (v) (2.12) 



2 



(h\ Mh R l / M*i R i 



3 



(h) MQ R M%M2 R 1 [MQ R , 

for the hierarchical scaling (this is the usual Edgeworth expansion) and 

p{ f \u,R) = ^H 3 (u) (2.13) 

p { /\v,R) = ^H 5 (u) 
p[ f \u,R) = I -!-,(— -) ) 




PS } <y,R) = ^ — 3^ — + ^tT) Hl{v 

for the feeder scaling. With this organization, the series Eq. (2.11) is an asymptotic expansion 
for PDFs with moments of either scaling. Organizing the expansion in this way also exposes 
the key difference between the hierarchical and feeder scalings: p^ 1 tx while p^' oc M^ 3 , 
which implies that higher-order terms in the expansion make a larger contribution to feeder 
distributions than they do to hierarchical distributions. 
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Then arranging the sums according to each scaling and performing the integral in Eq. 
(2.1), we find for objects of type X 



W> — (^) + ^ E E n £ <-> 

where the sets {fc m }/i are again non-negative integer solutions to fei + 2^2 + • ■ ■ + sk s = s and 

r = fei + ^2 + • • • + k m , but the sets {fc m } / are solutions to 3k± + 4^2 H h (s + 2)k s = s + 2. 

Now for either scaling, truncating the series at some finite s in the sums above keeps all terms 

s/3 

up to the same order in -M§ for hierarchical scalings and -M 3 for feeder scalings. 

We will work out constraints for the single parameter .M3 characterizing the two classes 
of distribution (hierarchical and feeder) we have motivated using their particular Petrov 
expansions shown above. First, we illustrate how the Petrov-type expansions work in Figure 
1, which shows a class of examples of the Edgeworth expansion compared to the corresponding 
complete distribution. The panels show three x 2 distributions with varying levels of non- 
Gaussianity and three truncations of the Edgeworth expansion for each case. As the figures 
show, the Edgeworth series is accurate out to larger values of v when the distribution is closer 
to Gaussian. But, for a given level of non-Gaussianity, there is a value of v beyond which 
the series is not an accurate fit and adding more terms does not improve the fit. A study 
comparing a variety of expansions of non- Gaussian distributions, with additional examples, 
can be found in [32]. 

2.2 Truncation and Error Evaluation 

To calculate the abundance of rare objects, we must truncate the f3 expansions given by 
Eq. (2.14), and we must determine the accuracy of the truncated series. For a given level of 
non-Gaussianity, there will be some region at high v beyond which the truncated expansion 
for /3 is not reliable. Furthermore, for a given value of v, there is also a level of non-Gaussianity 
beyond which the expansion is not reliable. This maximum accessible value of M3 can be 
increased by including higher-order terms, but only up to a point. Eventually the inclusion 
of more terms will stop improving the accuracy of the truncated sum, and adding additional 
terms makes the expansion less reliable (e.g., see the blue long-dashed curves in Figure 1). 
In this section, we determine the optimal number of terms to include in the expressions for 
the PDF and f3. We also describe how we evaluate the accuracy of the resulting expressions 
and find the maximum value of M3 that can be reliably constrained by a given upper bound 
on (3. Some previous work on the validity of truncating this series at low order and applying 
it to cluster number counts can be found in [33]. 

To calculate the number of objects of a given type, we must be able to trust the PDF 
for values of v > v Cj x that contribute to the f3 integral given by Eq. (2.1). Fortunately, 
only a limited range of v values make a significant contribution to (3 because the PDF 
decreases rapidly as v increases (~ exp[— u 2 /2\). Therefore, we only need to trust the PDF 
for v Cj x < v < z/ ma x) where f m ax is chosen so that truncating the j3 integral at ^ max has a 
minimal impact on (5. More explicitly, we demand that 

§p^4 < 0-02 (2.15) 

Pl^c.x) 
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Figure 1. A comparison of the Edge-worth expansion to the exact PDF for the xt distribution with 
mean zero and variance one, P(v) = V2k{V2kv + k) (k / 2 -^Exp[~(V2k v + fc)/2]/(2 fe / 2 r(fc/2)). The 
upper left panel shows the full distribution (solid black line) for k = 2, as well as the Edgeworth 
expansion truncated after p% (dotted red), p§ (short-dashed pink), and p-j (long-dashed blue), using 
hierarchical ordering. The upper right panel shows the same curves for k = 4, while the lower panels 
are for k — 20. The lower right is a zoomed-in view of the tail of the k — 20 distribution. 



so that the PDF for values of v > u max contributes less than two percent of the /3 integral. 
We find that this condition is satisfied for all interesting values of v C) x and M3 (i.e. f Cj x 
values for which f3 is large enough to be constrained and M.3 values for which we can trust 
the PDF for v Cj x < v < i/ max ) if we set ^ max = ifi^'}, where K = 2.1 for feeder models and 
K = 2.2 for hierarchical models. 
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Next, we evaluate the error in the PDF over the range v c x < v < f max as a function 
of M.3 and f Cj x • Since the series expansion of the PDF is an asymptotic expansion, the 
error introduced by truncating the series at N terms is the same order as the N + 1 term in 
the expansion. 4 Given M.3 and v c x-> we compute the ratio of the N + 1 term in the PDF 
to the sum of the first N terms and find its maximum value in the range v Ct x < v < f max : 
this is the maximum fractional error in the truncated PDF over the range of v values that 
contribute significantly to (3. In our analysis, we restrict ourselves to values of M3 such that 
this fractional error is less than twenty percent. Given that the PDF is accurate to within 
20% for all v values that significantly contribute to /3, we know that the fractional error in 
/3 is at most 20% for these M3 values. We expect that error in (5 is actually lower than 20% 
because the errors in the PDF are less than 20% for most values of v between v Ct x and f m ax- 
Furthermore, these errors may be positive or negative and therefore may partially cancel 
when the PDF is integrated to obtain /3. 

This expectation was confirmed when we applied our procedure to the xt probability 
distributions and compared the values of j3 obtained by integrating the truncated PDF ex- 
pansion over all v > u c to the values obtained by integrating the true PDF (given in the 
caption to Figure 1). We tested series that were truncated at s = 7,8,9, and 10 and we 
considered the series trustworthy only if the estimated PDF errors were less than 20% within 
the range v c x < v < Z'Wx- We estimated the errors in the PDF using the first term not 
included in the expansion; we did not compare the truncated PDF to the true PDF. For all 
truncations, the (3 values obtained from the truncated series with v c and k values in that 
truncation's trusted region differed from the true f3 values by less than 6%. 

Figure 2 shows the region of the u C; x~-M.3 plane that is excluded because the PDF 
errors are larger than 20% within the range v C) x <v < ^max for N = 12, 14, and 16 for both 
hierarchical and feeder models, and N = 17 for feeder models; the excluded regions lie above 
the thick curves. For hierarchical models, we see that adding terms with N > 12 does not 
significantly change the upper bound on computationally accessible M3 values for v c ^x < 8. 
For larger values of v Ct x, adding more terms to the PDF brings larger values of M3 into the 
trustworthy region. For feeder models, adding more terms decreases the range of accessible 
M3 values for small values of v c ,Xi but improves the range for larger values of v c ,x- However, 
for u c x > 8, adding more terms does not change the maximum acceptable value of M3. 

To determine the optimal number of terms to consider when computing the PDF, we 
must know what f Cj x values we are interested in. The shaded regions in Figure 2 show 
the regions that are excluded by different observational upper bounds on f3. To maximize 
the power of a given upper bound on j3, we want to choose ./V to maximize the area of the 
corresponding shaded region that lies underneath the computational constraint. For example, 
the blue and purple shaded areas show the values excluded by the bound f3 < 10~ 10 . For 
hierarchical models, we should choose N > 13 when evaluating this constraint. For feeder 
models N = 17 is the best choice because it maximizes the excluded range of u Cj x values, but 
we could reach higher M.3 by also considering PDFs with fewer terms. In practice, however, 
changing the number of terms in the expansion for j3 only changes the excluded regions for 
values of v C) x near the excluded region's boundary and for M.3 values that are far above all 
the thick curves in Figure 2. Therefore, we should choose N such that the computationally 

4 For both the hierarchical and the feeder models, we define the iVth term in the expansion according to 
the power of Ms in that term. For hierarchical models, the iVth term is proportional to and corresponds 
to the a — N term in Eq. (2.11). For feeder models, the A?th term is proportional to M^ 3 , and receives 
contributions from all terms with s + 2r = N in Eq. (2.11). 
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Figure 2. Computational and observational constraints for hierarchical models (left) and feeder 
models (right). The curves show the values of M3 above which the errors in the PDF truncated at 
N terms exceed 20% within the range of v values that contribute to /3(u Ct x)- The N for each curve is 
indicated directly on the curves. The shaded regions are excluded by observational constraints on /3 
with TV = 16 for hierarchical models and N = 17 for feeder models, assuming (3 < 10~ 6 , 10~ 10 , 10~ 15 
or 10~ 20 . Only the parts of the shaded regions falling below the colored curves are reliably excluded. 



accessible excluded region reaches the highest possible values of v c ,x- We can also trust this 
excluded region for all .M3 that are computationally accessible for any smaller value of N. 
Our default value for N is 16 for hierarchical models and 17 for feeder models. We chose 
this value of N due to computational limitations, but from Figure 2 we see that including 
more terms offers no benefits for hierarchical models with upper bounds on (3 greater than 
1CP 15 . For feeder models, including more terms could marginally extend the computationally 
accessible region excluded by (3 < 10 -10 , but it would not have any impact for lower or higher 
bounds on /3. 

Figure 3 compares several Gaussian and non- Gaussian PDFs (left panel) and their 
positive tails (right panel). The distributions have parameters found using the procedure 
above and the bound on the abundance of UCMHs at k = 2 x 10 7 Mpc _1 , for which 5 C ~ 
0.00156 and (3 < 3.0 x 10~ 6 . The thick solid black curve is the Gaussian distribution that 
saturates the bound on /3, with an = 0.000335. The blue solid curve is the non-Gaussian 
distribution with hierarchical scaling and the maximum trustworthy value of M3, -M3 imax 
0. 1 for iV = 12, that gives (3 = 3.0 x 10~ 6 . This case also defines a minimum variance 
or = 0.000259 for a hierarchical distribution that saturates the bound on j3; the dotted blue 
curve is the Gaussian distribution with that same value of a. In the right hand panel, this 
Gaussian distribution has very little area under the positive tail and lies nearly on top of the 
x-axis. The long-dashed red curve is the non-Gaussian distribution with feeder scaling and 
the maximum trustworthy value of .M3, -M3 jmax = 0.03 for N = 13, that gives (3 = 3.0 x 10~ 6 . 
This curve has <tr = 0.000308, and the dotted red curve in the right panel shows the tail of 
a Gaussian distribution with the same value of a. In the right panel, this curve would be 
nearly indistinguishable from the red dashed curve. The dot-dashed curve in the right panel 
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Figure 3. Comparing Gaussian and non- Gaussian PDFs for UCMHs at k = 2.0 x 10 7 Mpc _1 ; the 
right panel shows 6 > S c ~ 1.56x 10 -3 . In both panels, the thick solid black curve is the Gaussian PDF 
that saturates the observational bound on the UCMH abundance (j3 < 3.0 x 10 -6 ). Both panels also 
show the maximally non-Gaussian hierarchical PDF (solid blue) and feeder PDF (dashed red) that 
saturate this bound. (That is, with the maximum value of M3 for which the expansion is controlled 
out to the largest value of v that contributes significantly to the integral.) In the right panel, we 
indicate that all three of these curves have the same value of (3. The dotted curves show Gaussian 
PDFs that have the same values of a as these maximally non-Gaussian PDFs (blue for hierarchical 
and red for feeder). In addition, the dot-dashed curve in the right hand panel shows a hierarchical 
PDF with the same values of a and M3 as the maximally non-Gaussian feeder PDF that saturates 
the (3 bound (the red dashed curve). 



is one way of visualizing the difference between the feeder and hierarchical scalings: it shows 
a hierarchical non-Gaussian PDF with values of a and .M3 that saturate j3 when used in a 
feeder PDF (M.3 )Xna x = 0-03 and or = 0.000308, as shown by the red dashed curve). Clearly, 
for a particular value of .M3, the feeder scaling gives a PDF that is overall substantially more 
non-Gaussian than the hierarchical scaling. 

Now we can get some measure of the maximum shift in the PDF that can come from 
adding weak non-Gaussianity in a controlled way. Figure 4 illustrates this in two ways. The 
left hand panel illustrates how a constraint on the variance o derived assuming a Gaussian 
distribution can be shifted to a constraint on a non-Gaussian distribution with a new (smaller) 
variance. The vertical axis shows the relative shift in v c between the Gaussian distribution 
that gives abundance fraction fj and the maximally non-Gaussian distribution, under control 
according to our criteria above, that gives the same fi. To use this for any particular object, 
v c should be written in terms of the appropriate S c and or for the object of interest. -M3 !max 
can be likewise be converted to a non-Gaussian parameter in the primordial distribution. 
This plot gives a sense of how accurately the variance of fluctuations must be known in order 
for number counts to provide an independent constraint on non-Gaussianity using Petrov 
expansions. For instance, if Av/vq = 0.3 in the left panel of Fig. 4, then adding weak 
non-Gaussianity in a controlled way has the same effect on the abundance of rare objects as 
increasing or by ~30%. If the uncertainty in or exceeds 30%, then number counts with the 
corresponding value /3 cannot independently constrain non-Gaussianity without knowing the 
full PDF of the density fluctuations. 

The right hand panel of Figure 4 shows for a given object and fixed variance (so fixed u c ), 
by what maximum factor the abundance f3 can change by adding controlled non-Gaussianity, 
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Figure 4. Left panel: The ratio [y c ,NG(-^3,max) — ^c.gI/^Gj where j^g is the value of v c such that 
a Gaussian distribution gives the plotted value of (3 (see Eq. (2.2)) and i / c ,NG(-^ / '3,max) is the value of 
v c such that a non-Gaussian distribution with Ai^ = A^max gives the same value of ft. Here A^max 
is the largest amplitude of non-Gaussianity where the Petrov expansion, calculated out to N = 16 for 
hierarchical scaling and N = 17 for feeder scaling, is trustworthy, according to the criteria discussed 
below Eq. (2.15). Right panel: The ratio of [/3(Ai3 jlna x, v c ) — ^c)]//3(0, v c ) plotted versus v c 
assuming hierarchical (blue circles) or feeder (red squares) scaling. In both plots, the change in slope 
in the hierarchical points appears to be an artifact of our truncation to a finite number of terms - going 
beyond N — 16 would probably allow one to continue go farther onto the tail. For the feeder scaling, 
however, it is less clear that adding more terms will help - the distribution is quite non-Gaussian, so 
the Petrov expansion breaks down irreparably at smaller v and smaller A13. The scatter in the points 
about a smooth line comes from the approximate nature of our criteria for determining A^max- 



[/3(A43 jmax , u c )—(3(0, u c )]/ /3(0, u c )- Notice that the primary limit on how much f3 can change is 
our ability to accurately determine the PDF. Since the feeder scaling produces a distribution 
that is overall much more non-Gaussian than the hierarchical for a fixed value of A^3, we 
cannot accurately calculate at as large of v c (rarer objects) or for as large of shifts in f3 as for 
the hierarchical case. Again, the effects of more non-Gaussian distributions can be computed 
if one knows the full distribution rather than just a finite number of moments. 

3 Primordial Black Holes 

For a Gaussian spectrum of fluctuations, the present-day fractional density of black holes has 
typically been used to constrain the spectral index n s under the assumption of a scale-free 
spectrum [14, 18, 34, 35]. The resulting bound, n s < 1.3 (based on e.g. VERITAS searches 
at k = 3 x 10 15 Mpc -1 [36]), is easily satisfied by the current experimentally measured 
value n s = 0.968 ± 0.012 [1], with little or no running. In terms of generalized, scale- 
dependent Gaussian power spectra, PBHs have been used to place bounds on the amplitude 
of curvature perturbations at the level of V-ji < 10~ 2 [8, 37], on scales corresponding to 
10~ 2 < k < 10 19 Mpc~ 1 . Stronger constraints than this can be found if the analysis is done 
using a top-hat windowing function rather than a Gaussian, which is arguably a more correct 
method [11, 38]. 

Present observational bounds on PBHs [8, 9] limit Ppbh to values ranging from fipBH ^5 
10~ 10 for solar mass PBHs, to (3pbh ^ 10 -28 for a small range of PBH masses around 10 13 g, 
to Ppbh ^ 10 -20 for even smaller masses. A PBH mass Mpbh corresponds to horizon mass 
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Mh = 3 3//2 Mpbh [8], and therefore wavenumber 



k 



eq 




(3.1) 



Here k eq = 0.072O m /i 2 Mpc" 1 = 9.68 x lO^Mpc^ 1 [1] and Af£ q = 3.5 x 10 17 M [13] are 
the wavenumber and horizon mass respectively at equality, and g 1 ^ = 3.91, gl q = 3.36, 
> 106.75 [39] are the respective degrees of freedom at equality and the time of horizon 
entry. 

We use the current limits on (3pbh from [9], along with Eq. (2.14) and the requirement 
that it represents a well-behaved expansion, to determine the range of values of Ms(k) that 
are excluded by non-observation of PBHs. For this calculation we draw upon the expressions 
for or in a generalized Gaussian power spectrum, given in Appendix B of [11]; in particular, 
we employ the top- hat window function advocated in that paper. We give our limits in Fig. 5 
as a function of the power spectrum of curvature perturbations at any given k 



We show the excluded areas as dark shaded regions, for both hierarchical and feeder scaling, 
and at three different example scales. The corresponding limits on (3 at these scales are 
P < 3.0 x 10~ 12 (k = 2 x 10 3 Mpc~ 1 ), < 1.9 X 10~ n (k = 2 x 10 7 Mpc~ 1 ) and /3 < 
3.5 x 10~ 28 (k = 2 x 10 17 Mpc _1 ). More positive non-Gaussianity (larger M.%) results in 
stronger exclusions, up until the point at which our expansion Eq. (2.14) begins to break 
down. Above these very large values of .M3, the excluded regions shrink with additional non- 
Gaussianity, as more of the available parameter space is made inaccessible by the increasingly 
poor behavior of the expansion. 

4 Ultracompact Minihalos 

UCMHs were initially discussed in early papers by Berezinsky and collaborators [40-43], 
and have been recently studied in the context of lensing signatures [12, 27], gamma-ray and 
neutrino signals from dark matter annihilation or decay [10, 11, 13, 44-48] and impacts on 
the CMB [48-51]. To date the strongest limits come from gamma-ray observations [11], 
indicating that n s < 1.17 for a scale-free spectrum, and V-n(k) < 10 _6 -10 -7 for generalized 
spectra, on scales corresponding to 3Mpc _1 < k < 3 x 10 7 Mpc _1 . The corresponding direct 
limits on j3 at the two examples scales of Fig. 5 accessible with UCMHs are f3 < 1.8 x 10~ 9 
(Jfc = 2 x lO^pc" 1 ) and < 3.0 x 10~ 6 (k = 2 x 10 7 Mp C - 1 ). 

To explore the impact of UCMHs on the allowed parameter space for non-Gaussian- 
ity, we have implemented Eq. (2.14) in the analysis framework of [11], using combined 1- 
year Fermi-LAT limits on the abundance of UCMHs from Galactic and extragalactic source 
searches 5 , as well as the overall diffuse gamma-ray flux in the direction of the Galactic poles. 
Whilst 2-year source searches (e.g. [52]) and corresponding diffuse data are available, we 
conservatively adopt 1-year sensitivities to be certain that no dark matter sources were found 
in the observing period we consider. The reader is referred to [11] for extensive discussion of 
this issue. 



5 See Appendix A for details of our slight improvement of the treatment of sources over [11]. 



Pn(k) = 1.104. 



(3.2) 
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Figure 5. 95% CL bounds on M.%(k) from UCMHs and PBHs, as a function of the power spectrum 
of curvature perturbations. The limits are given at k = 2 x 10 3 Mpc _1 (top), k = 2 x 10'Mpc -1 
(center) and k = 2 x 10 17 Mpc _1 (bottom), for both the hierarchical and feeder scaling scenarios. In 
general UCMHs probe non-Gaussianity at much lower overall amplitudes, but PBHs are sensitive to 
a much broader range of scales. 
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Our gamma-ray limits assume the canonical thermal cross-section for dark matter anni- 
hilation (3 x 10 -26 cm 3 s _1 ) into b quarks and anti-quarks, and a relatively conservative dark 
matter mass of 1 TeV (smaller masses would result in stronger limits). In Fig. 5, we show the 
resulting excluded values of Ai^ for the two scaling scenarios as light shaded regions, under 
different assumptions about the power spectrum of perturbations. As expected, the UCMH 
limits are considerably stronger than those from PBHs and show the same contraction at 
extremely large .M3, due to the breakdown of the expansion Eq. (2.14) in this regime. They 
do not provide any constraint on the smallest scales, however. They also do not reach the 
level of power expected from a simple scale-invariant extrapolation of the spectrum observed 
on CMB scales, which is V n ~ 2 x 1(T 9 [2]. 

The 95% confidence level (CL) exclusion regions we have plotted in Fig. 5 assume that 
only objects undergoing gravitational collapse before a redshift of z c = 1000 form UCMHs. 
This follows earlier work [11, 27] and is a very conservative choice. The defining characteristic 
of a UCMH is that it forms in isolation from a spherical collapse; the processes of radial col- 
lapse and secondary infall give UCMHs extremely steep dark matter density profiles, thereby 
distinguishing them from halos formed at late times. Exactly when the approximation of 
radial infall breaks down is not yet known [11], as this depends sensitively on the velocities 
of dark matter particles and interactions between growing overdensities; more detailed simu- 
lation work is sorely needed to answer this question. In principle UCMHs may form at lower 
redshifts, as long as the presence of different-sized overdensities do not spoil the assumption 
of spherical collapse; the smaller the redshift, the more unlikely it is that this requirement 
will be fulfilled. In Fig. 5 we plot the excluded values of .A/f 3 as a function of Vn, under the 
assumption that z c = 200. Later collapse redshifts allow smaller amplitude perturbations to 
form UCMHs, reducing <5 C) ucmh and increasing /?ucmh for a given combination of M.3 and 
Viz- If UCMHs continue to be formed as late as z c = 200, their non-observation begins to 
constrain non-Gaussianities on small scales, even if the power spectrum is near the equivalent 
scale-invariant level. 

5 Discussion &; Conclusions 

We have provided a general technique for using the abundance of any object to constrain 
weakly non-Gaussian statistics of the primordial fluctuations. We suggest precise criteria 
for the amplitude of non-Gaussianity that can be accurately constrained using Petrov-type 
expansions and emphasize that there is not a one-to-one mapping between the overall level 
of non-Gaussianity and the size of the skewness (.M3). This is important when using the 
abundance of rare objects to constrain non-Gaussianity: number counts are sensitive to any 
deviation from Gaussianity, but a measured value of (3 does not constrain a parameter like 
Inl (typically interpreted as the amount of skewness) without additional information about 
the higher moments of the PDF. 

In general PBHs provide stringent constraints on non-Gaussianities only in the presence 
of fluctuations with relatively large variance, from V-ji > 1 - 2 x 10~ 3 at intermediate {k = 
2 x 10 3 Mpc -1 ) and small (k = 2 x 10 7 Mpc -1 ) scales, to Vn >7x 10~ 4 at very small scales. 
Notably though, at very small scales PBHs provide the only available limit on the allowed 
amount of non-Gaussianity. UCMHs provide stronger constraints, and so may be more likely 
to eventually also constrain non-Gaussianity, but do not reach as small scales as PBHs. 

For both UCMHs and PBHs, it is straightforward to map our constraints on M3 to 
constraints on a non-Gaussian parameter in the primordial fluctuations. For example, for 
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the local ansatz lZ(x) = lZ(x)c + \fNh[R-G{x) 2 — (R-g(x) 2 }] the smoothed, integrated, dimen- 
sionless skewness A4% ~ (1Z(x) 3 } r/ \{1Z{x) 2 ) i?,) 3 ^ 2 is very nearly equal to /nl'PjI 2 - For this 
case, if the variance is large enough, the abundance of UCMHs constrains Jnl ^ 0(100 ) and 
PBHs constrain /jvl < O(10) on smaller scales. However, it is important to remember that 
this does not mean UCMHs are less powerful - quite the opposite. For the constraint from 
either object, the maximum value of Jnl corresponds to M.3 ~ O(0. 1), so those values of 
Inl correspond to equally non-Gaussian distributions. The difference in constraints phrased 
as a value of /jvx comes from the difference in the variance of the fluctuations. For refer- 
ence, recall that on CMB scales, where the amplitude of fluctuations is currently measured, 
f NL ~ O(10) corresponds to Ma ~ O(10~ 3 ). 

To constrain non-Gaussianity, both PBHs and UCMHs require significantly more power 
than extrapolating the CMB spectrum would suggest, but PBHs require even more power 
than UCMHs. So, given a prediction for the power on small scales, UCMHs are more likely to 
constrain non-Gaussianity. However, these constraints assume that the dark matter particle 
is a standard thermal relic with detectable annihilation products [11]. Unfortunately, near- 
future astrometric microlensing searches for UCMHs can only constrain /3ucmh ~ 0.3 [12]. 
For these large values of /3, UCMHs are not sufficiently rare to be a useful probe of non- 
Gaussianity because v c ~ 1, which is smaller than the crossing point between Gaussian and 
skewed PDFs. Consequently, changing M3 does not significantly affect the abundance of 
relatively common objects (f3 > 0.1). 

Aside from the abundances of rare objects, the only other probe of the small-scale 
inhomogeneities is the spectral distortion of the CMB, specifically the //-distortion, which is 
sensitive to the fluctuations on scales on the order of k = 50 — 10 4 Mpc -1 [15-17]. As with 
PBHs and UCMHs, CMB //-distortion can also constrain non-Gaussianity, but the nature of 
the constraints on both characteristics is qualitatively different. First, CMB //-distortion is 
sensitive to the total power in fluctuations over a wider range of scales, while the objects are 
sensitive to the power only very near the scale of the object, k ~ 1/R. Second, //-distortion is 
sensitive to only the squeezed limit of the bispectrum rather than the total skewness. That is, 
this measurement can constrain the amplitude in the correlation of three momentum modes 
with ki -C &2 ~ ks [53, 54]. If non-Gaussianity follows the local ansatz (which has a strong 
signal in the squeezed bispectrum), a futuristic CMB probe like PIXIE [55] could constrain 
Inl ~ O(10 3 ) for a scale-invariant power spectrum and /jvl ~ O(10) for an enhanced 
small-scale power spectrum that saturates current bounds (M3 ~ 0.04 in both cases) [53]. 
This probe is then quite complementary to the abundance of PBHs and UCMHs in terms of 
the range of scales probed and the relationship between the observables and the primordial 
fluctuations. 
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A Generalized source limits on UCMHs 



Our analysis here actually includes a slightly improved form of the point source analysis of 
[11]. Following from Eq. (18) in [11], if there is a single UCMH present somewhere in the 
Milky Way, we identify the probability P 1 seen of actually seeing that UCMH as 

pseen = p ^exists = jojoba r A y 

Here d is the true distance of the UCMH from Earth, (i b s is the maximum observable 
distance required for detection at CL x, P^^f i is the probability of the UCMH residing 
within distance (i b s of Earth, Mmw is the dark mass of the Milky Way and M^ <( ^ obs is the 
dark mass within d bs- 

This leads to an improved version of Eq. (24) in [11], which gives the maximum allowed 
fraction of matter in UCMHs today 



f x<CMK log(l-l/) 



f _ JX J "UCMH ^6V^ yj /a 9 n 



Here f x = p x / p m is the fraction of matter in the Universe that is dark, Af{j CMH is the present- 
day UCMH mass, and y is the desired CL of the limit / max - The relationship between j3 and 
/max can be found in Eq. (21) of [11]. Eq. (A. 2) is valid for all x,y < 1. For 1 — y S> 1 — x, 
it agrees with Eq. (24) in [11]; whilst both forms are valid for the values of x and y we use 
here (1 — y = 5 x 1(T 2 , 1 - x = 6 x 1CT 7 , as in [11]), the latter breaks down as y — > x, so Eq. 
(A. 2) is more general. 



References 



[1] E. Komatsu et. at, Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 
Cosmological Interpretation, Astrophys.J.Suppl. 192 (2011) 18, [arXiv: 1001 .4538]. 

[2] J. Dunkley, R. Hlozek, J. Sievers, V. Acquaviva, P. Ade, et. al, The Atacama Cosmology 

Telescope: Cosmological Parameters from the 2008 Power Spectra, Astrophys.J. 739 (2011) 52, 
[arXiv: 1009.0866]. 

[3] R. Keisler, C. Reichardt, K. Aird, B. Benson, L. Bleem, et. ai, A Measurement of the Damping 
Tail of the Cosmic Microwave Background Power Spectrum with the South Pole Telescope, 
Astrophys.J. 743 (2011) 28, [arXiv: 1105.3182]. 

[4] B. A. Reid, W. J. Percival, D. J. Eisenstein, L. Verde, D. N. Spergel, et. at, Cosmological 
Constraints from the Clustering of the Sloan Digital Sky Survey DR7 Luminous Red Galaxies, 
Mon.Not.Roy.Astron.Soc. 404 (2010) 60-85, [arXiv : 0907 . 1659]. 

[5] N. Sehgal, H. Trac, V. Acquaviva, P. A. Ade, P. Aguirre, et. at, The Atacama Cosmology 
Telescope: Cosmology from Galaxy Clusters Detected via the Sunyaev-Zel'dovich Effect, 
Astrophys.J. 732 (2011) 44, [arXiv: 1010. 1025]. 

[6] SDSS Collaboration Collaboration, P. McDonald et. ai, The Lyman- alpha forest power 
spectrum from the Sloan Digital Sky Survey, Astrophys.J.Suppl. 163 (2006) 80-109, 
[astro-ph/0405013] . 



- 16 - 



[7] B. J. Carr, The primordial black hole mass spectrum, Astrophys. J. 201 (1975) 1-19. 

[8] A. S. Josan, A. M. Green, and K. A. Malik, Generalised constraints on the curvature 

perturbation from primordial black holes, Phys.Rev. D79 (2009) 103520, [arXiv: 0903. 3184]. 

[9] B. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama, New cosmological constraints on primordial 
black holes, Phys.Rev. D81 (2010) 104019, [arXiv: 0912. 5297]. 

[10] A. S. Josan and A. M. Green, Gamma-rays from ultracompact minihalos: potential constraints 
on the primordial curvature perturbation, Phys.Rev. D82 (2010) 083527, [arXiv: 1006.4970]. 

[11] T. Bringmann, P. Scott, and Y. Akrami, Improved constraints on the primordial power 
spectrum at small scales from ultracompact minihalos, Phys.Rev. D85 (2012) 125027, 
[arXiv: 1110.2484]. 

[12] F. Li, A. L. Erickcck, and N. M. Law, A new probe of the small-scale primordial power 

spectrum: astrometric microlensing by ultracompact minihalos, Phys.Rev. D86 (2012) 043519, 
[arXiv: 1202. 1284]. 

[13] P. Scott and S. Sivertsson, Gamma-Rays from Ultracompact Primordial Dark Matter 
Minihalos, Phys.Rev. Lett. 103 (2009) 211301, [arXiv: 0908. 4082]. 

[14] J. S. Bullock and J. R. Primack, Non-gaussian fluctuations and primordial black holes from 
inflation, Phys. Rev. D55 (1997) 7423-7439, [astro-ph/9611106]. 

[15] W. Hu, D. Scott, and J. Silk, Power spectrum constraints from spectral distortions in the 
cosmic microwave background, Astrophys. J. 430 (1994) L5-L8, [astro-ph/9402045]. 

[16] J. Chluba and R. A. Sunyaev, The evolution of CMB spectral distortions in the early Universe, 
MNRAS 419 (2012) 1294-1314, [arXiv: 1109.6552]. 

[17] J. Chluba, A. L. Erickcek, and I. Ben-Dayan, Probing the inflaton: Small-scale power spectrum 
constraints from measurements of the CMB energy spectrum, Astrophys. J. 758 (2012) 76, 
[arXiv: 1203.2681]. 

[18] A. M. Green and A. R. Liddle, Constraints on the density perturbation spectrum from 
primordial black holes, Phys. Rev. D56 (1997) 6166-6174, [astro-ph/9704251]. 

[19] P. Ivanov, Nonlinear metric perturbations and production of primordial black holes, Phys.Rev. 
D57 (1998) 7145-7154, [astro-ph/9708224]. 

[20] P. Pina Avelino, Primordial black hole constraints on non-gaussian inflation models, Phys.Rev. 
D72 (2005) 124004, [astro-ph/0510052]. 

[21] J. C. Hidalgo, The effect of non- Gaussian curvature perturbations on the formation of 
primordial black holes, arXiv : 0708 . 3875. 

[22] P. Klimai and E. Bugaev, Primordial black hole formation from non-Gaussian curvature 
perturbations, arXiv: 1210.3262. 

[23] C. T. Byrnes, E. J. Copeland, and A. M. Green, Primordial black holes as a tool for 
constraining non-Gaussianity, Phys.Rev. D86 (2012) 043512, [arXiv: 1206.4188]. 

[24] D. S. Salopek and J. R. Bond, Nonlinear evolution of long wavelength metric fluctuations in 
inflationary models, Phys. Rev. D42 (1990) 3936-3962. 

[25] W. H. Press and P. Schechter, Formation of galaxies and clusters of galaxies by self similar 
gravitational condensation, Astrophys. J. 187 (1974) 425-438. 

[26] M. LoVerde, A. Miller, S. Shandera, and L. Verde, Effects of Scale-Dependent Non-Gaussianity 
on Cosmological Structures, JCAP 0804 (2008) 014, [arXiv: 0711 . 4126]. 

[27] M. Ricotti and A. Gould, A New Probe of Dark Matter and High-Energy Universe Using 
Microlensing, Astrophys. J. 707 (2009) 979-987, [arXiv: 0908. 0735]. 



-17- 



[28] V. V. Petrov, Sums of independent random variables. Springer- Verlag, Berlin, New York, 1975. 

[29] V. V. Petrov, Limit theorems of probability theory: sequences of independent random variables, 
vol. 4 of Oxford studies in probability. Clarendon Press, Oxford / Oxford University Press, New 
York, 1995. 

[30] N. Barnaby and S. Shandera, Feeding your Inflaton: Non-Gaussian Signatures of Interaction 
Structure, J CAP 1201 (2012) 034, [arXiv : 1109 . 2985]. 

[31] X. Chen and Y. Wang, Quasi-Single Field Inflation and Non-Gaussianities, JCAP 1004 (2010) 
027, [arXiv: 091 1.3380]. 

[32] S. Blinnikov and R. Moessner, Expansions for nearly Gaussian distributions, 
Astron.Astrophys.Suppl.Ser. 130 (1998) 193-205, [astro-ph/9711239]. 

[33] S. Chongchitnan and J. Silk, A Study of High- Order Non-Gaussianity with Applications to 
Massive Clusters and Large Voids, Astrophys.J. 724 (2010) 285-295, [arXiv : 1007 . 1230]. 

[34] A. M. Green, A. R. Liddle, K. A. Malik, and M. Sasaki, A new calculation of the mass fraction 
of primordial black holes, Phys. Rev. D70 (2004) 041502, [astro-ph/040318l]. 

[35] S. Chongchitnan and G. Efstathiou, Accuracy of slow-roll formulae for inflationary 

perturbations: Implications for primordial black hole formation, JCAP 0701 (2007) 011, 
[astro-ph/0611818]. 

[36] G. Tesic, for the VERITAS Collaboration, Searching for primordial black holes with the 
VERITAS gamma-ray experiment, J.Phys.Conf.Ser. 375 (2012) 052024. 

[37] L. Alabidi, K. Kohri, M. Sasaki, and Y. Sendouda, Observable Spectra of Induced Gravitational 
Waves from Inflation, JCAP 1209 (2012) 017, [arXiv: 1203.4663]. 

[38] T. Bringmann, C. Kiefer, and D. Polarski, Primordial black holes from inflationary models with 
and without broken scale invariance, Phys. Rev. D65 (2002) 024008, [astro-ph/0109404]. 

[39] E. W. Kolb and M. S. Turner, The Early Universe. Frontiers in Physics, Addison- Wesley, 
Reading, MA, 1990. 

[40] V. Berezinsky, V. Dokuchaev, and Y. Eroshenko, Small - scale clumps in the galactic halo and 
dark matter annihilation, Phys. Rev. D68 (2003) 103003, [astro-ph/0301551]. 

[41] V. Berezinsky, V. Dokuchaev, and Y. Eroshenko, Destruction of small-scale dark matter clumps 
in the hierarchical structures and galaxies, Phys. Rev. D73 (2006) 063504, [astro-ph/0511494]. 

[42] V. Berezinsky, V. Dokuchaev, and Y. Eroshenko, Anisotropy of dark matter annihilation with 
respect to the Galactic plane, JCAP 0707 (2007) 011, [astro-ph/0612733]. 

[43] V. Berezinsky, V. Dokuchaev, and Y. Eroshenko, Remnants of dark matter clumps, Phys. Rev. 
D77 (2008) 083519, [arXiv: 0712. 3499]. 

[44] B. C. Lacki and J. F. Beacom, Primordial Black Holes as Dark Matter: Almost All or Almost 
Nothing, Astrophys.J. 720 (2010) L67-L71, [arXiv: 1003.3466]. 

[45] V. Berezinsky, V. Dokuchaev, Y. Eroshenko, M. Kachclriess, and M. A. Solberg, Superdense 
cosmological dark matter clumps, Phys. Rev. D81 (2010) 103529, [arXiv : 1002 . 3444]. 

[46] V. Berezinsky, V. Dokuchaev, Y. Eroshenko, M. Kachelriess, and M. A. Solberg, Annihilations 
of superheavy dark matter in superdense clumps, Phys. Rev. D81 (2010) 103530, 
[arXiv: 1002.3445]. 

[47] Y.-P. Yang, L. Feng, X.-Y. Huang, X. Chen, T. Lu, et. al, Constraints on ultracompact 

minihalos from extragalactic j-ray background, JCAP 1112 (2011) 020, [arXiv : 1112 . 6229]. 

[48] Y. Yang, G. Yang, and H. Zong, Dark Matter Decay and the Abundance of Ultracompact 
Minihalos, arXiv: 1210. 1409. 



- 18 - 



[49] D. Zhang, Impact of primordial ultracompact minihaloes on the intergalactic medium and first 
structure formation, MNRAS 418 (2011) 1850-1872, [arXiv: 1011 . 1935]. 

[50] Y. Yang, X. Huang, X. Chen, and H. Zong, New Constraints on Primordial Minihalo 

Abundance Using Cosmic Microwave Background Observations, Phys.Rev. D84 (2011) 043506, 
[arXiv: 1109.0156]. 

[51] Y.-P. Yang, X. Chen, T. Lu, and H.-S. Zong, The Abundance of New Kind of Dark Matter 
Structures, Eur.Phys.J .Plus 126 (2011) 123, [arXiv : 1112 . 6228]. 

[52] H.-S. Zechlin and D. Horns, Unidentified sources in the Fermi-LAT second source catalog: the 
case for DM subhalos, arXiv: 1210.3852. 

[53] E. Pajer and M. Zaldarriaga, A New Window on Primordial non-Gaussianity, Phys.Rev. Lett. 
109 (2012) 021302, [arXiv : 1201 . 5375]. 

[54] J. Ganc and E. Komatsu, Scale- dependent bias of galaxies and mu-type distortion of the cosmic 
microwave background spectrum from single-field inflation with a modified initial state, 
Phys.Rev. D86 (2012) 023518, [arXiv: 1204.4241]. 

[55] A. Kogut, D. Fixsen, D. Chuss, J. Dotson, E. Dwek, et. al., The Primordial Inflation Explorer 
(PIXIE): A Nulling Polarimeter for Cosmic Microwave Background Observations, JCAP 1107 
(2011) 025, [arXiv: 1105.2044]. 



- 19 - 



